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ABSTRACT 

In the core-accretion model, gas-giant planets form solid cores which then ac- 
crete gaseous envelopes. Tidal interactions with disk gas cause a core to undergo 
inward type-I migration in 10 4 to 10 5 years. Cores must form faster than this 
to survive. Giant planets clear a gap in the disk and undergo inward type-II 
migration in < 10 6 years if observed disk accretion rates apply to the disk as a 
whole. Type-II migration times exceed typical disk lifetimes if viscous accretion 
occurs mainly in the surface layers of disks. Low turbulent viscosities near the 
midplane may allow planetesimals to form by coagulation of dust grains. The 
radius r of such planetesimals is unknown. If r < 0.5 km, the core formation time 
is shorter than the type-I migration timescale and cores will survive. Migration is 
substantial in most cases, leading to a wide range of planetary orbits, consistent 
with the observed variety of extrasolar systems. When r ~ 100m and midplane 
a ~ 3 x 10 -5 , giant planets similar to those in the Solar System can form. 

Subject headings: planetary systems: formation — planetary systems: protoplan- 
etary disks — planets and satellites: formation — solar system: formation 



1. Introduction 

In the core-accretion model, gas-giant planets begin life as solid cores that grow by 
sweeping up small planetesimals (Inaba et al. 2003). Cores > 10~ 5 Af e grow oligarchically: 
each radial zone in a protoplanetary disk contains a single core and many planetesimals 
(Kokubo & Ida 1998; Thommes et al. 2003). Perturbations from a core control the orbital 
distribution of nearby planetesimals which then determines the core's growth rate. Cores 
larger than Mars acquire extended atmospheres of nebular gas. Passing planetesimals are 
slowed by gas drag as a result, increasing the chance of capture (Inaba & Ikoma 2003). When 
a core reaches a critical mass M crit ~ 10M e , a static atmosphere can no longer be supported, 
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and the core accretes a massive gas envelope, becoming a gas-giant planet (Hubickyj et al. 
2005; Ikoma et al. 2000). 

Cores must grow to M crit before the disk gas disperses, which typically happens in a few 
Myr (Haisch et al. 2001). A core of mass M sweeps up planetesimals at a rate given by 



where P and r h are the core's orbital period and Hill radius, and S so iid is the local surface 
density of planetesimals (Inaba et al. 2001). The collision probability P co i depends on the 
relative velocity v TC \ of passing planetesimals, which is a function of their mean eccentricity 
e and inclination % (Inaba et al. 2001). 

In the minimum-mass solar nebula (MMSN), the solid and gas surface densities are: 



where E r/i ce ~ 3 gem" 2 and x ~ —3/2 (Weidenschilling 1977). Here, a is the orbital radius; 
S is the surface density of rocky material at 5 AU; r/ gas ~ 200; rj lcc = 1 inside the snow 
line, and r] icc > 1 outside the snow line. Most models for the formation of Jupiter use 
mass-enhanced nebulae, with S 77 ice = 8-25 gem -2 , in order to grow to M crit before the gas 
disperses (Inaba et al. 2003; Alibert et al. 2005; Hubickyj et al. 2005; Thommes et al. 2003). 

A core generates spiral density waves in the gas which cause the core to undergo inward 
type-I migration at a rate given by: 



where v\ ex> is the core's Keplerian orbital velocity, c s is the gas sound speed, and M* is the 
stellar mass (Tanaka et al. 2002). (In this paper we consider only migration caused by the 
gaseous component of the disk). A 10M e core at 5 AU in a mass-enhanced nebula will 
migrate into its star in 10,000-30,000 years, much less than the time required to form a core 
in most models (Inaba et al. 2003; Alibert et al. 2005; Hubickyj et al. 2005; Thommes et al. 
2003; Ida & Lin 2004). Growth and migration rates are both proportional to E (for a given 
disk metallicity), so this result is independent of the disk mass. 

A massive planet clears an annular gap in the disk and undergoes type-II migration, 
moving inwards at the same rate that gas flows viscously towards the star: 
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where a = vv^j \ac 2 s ) and v is the disk viscosity (D'Angelo et al. 2002). Observed disk 
accretion rates imply a ~ 0.001 (Hueso & Guillot 2005), so a planet at 5 AU will migrate 
into its star in 0.5 Myr. This is less than the lifetime of most disks (Haisch et al. 2001). 

2. Living with Migration 

The existence of gas-giant planets suggests that: (i) giant-planet cores grow rapidly 
before type-I migration moves them into the star, and (ii) giant planets form at locations 
where the disk viscosity is lower than observed disk accretion rates would suggest. 

Core growth should be fastest when planetesimals are small for two reasons. Small plan- 
etesimals experience strong gas drag, reducing e and i sufficiently for v rc \ to be determined by 
Keplerian shear in the disk (Rafikov 2004; Chambers 2006). As a result, P co i is much higher 
than in the dispersion-dominated regime considered by most models. In addition, small 
planetesimals are slowed when they pass through the atmospheres of large cores, increasing 
the capture probability (Inaba & Ikoma 2003). 

Most previous models have considered planetesimals with radii r = 10-100 km (Inaba 
et al. 2003; Alibert et al. 2005; Hubickyj et al. 2005; Thommes et al. 2003). Planetesimals 
would have this size if they formed via gravitational instabilities (GI) in the solid component 
of the disk (Wetherill 1980). However, GI requires S so i id /E gas ~ 1 (Garaud & Lin 2004), 
which probably occurs under only limited circumstances. Alternatively, planetesimals may 
form by pairwise coagulation of dust grains (Weidenschilling 1997), but it is unclear how 
large these objects will be. Regardless of the initial planetesimal size, r probably becomes 
small once cores grow larger than Ceres, due to collisional fragmentation (Kenyon & Bromley 
2004). Given these uncertainties, we will treat r as a model parameter. 

We can gauge the importance of r using a simple model for core growth based on 
equation 1. Here, e and i are determined assuming an equilibrium between excitation due to 
perturbations from nearby cores (Ohtsuki et al. 2002) and damping due to gas drag (Inaba et 
al. 2001). The core's capture cross section is calculated including the effect of its atmosphere 
using equations A7-A12 of Inaba & Ikoma (2003). Core growth rates are assumed to be 50% 
higher than in equation 1 due to collisions between neighbouring cores (Chambers 2006). 

Figure 1 shows the mass of a core growing at 5 AU, neglecting migration, for 4 values 
of r. Here, £ S oiid = 10 gem -2 , solids have a density of 1.5 gem -3 , and E gas /S so iid = 90. Core 
growth times are shorter when r is small, as expected. When r < 100 m, the core formation 
time is shorter than the type-I migration timescale, shown by the dashed line. Dissipation 
of the nebular gas will slow type-I migration, and migration will speed up core growth by 
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increasing the supply of planetesimals (Alibert et al. 2005). Hence, cores should survive 
when r is somewhat larger than 100 m. 

The source of viscosity in protoplanetary disks is unclear. One possibility is magneto- 
rotational instability (Hawley et al. 1995), which is highly effective but probably confined 
mainly to the surface layers of a disk (Matsumura & Pudritz 2006). Near the midplane, where 
planets form, v ~ 30 x lower than in the surface layers (Turner et al. 2006), corresponding 
to a ~ 3 x 10^ 5 . Similar values of a can explain the size distribution of chondrules seen in 
meteorites (Cuzzi et al. 2001). Low a also promotes dust coagulation by reducing typical 
collision speeds and fragmentation. For a = 3x 10 -5 , the type-II migration time for a planet 
at 5 AU is ~ 15 Myr, longer than the lifetime of most disks (Haisch et al. 2001). 

A core opens a gap in the disk when M > M gap = max(M nv , M vis ), where M nv is the 
mass needed to open a gap in a non-viscous disk, given by 
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where Q = c s M+ / 'fk e pS gas 7ra 2 is the Toomre stability criterion (Rafikov 2002); and M v 
the mass needed to maintain a gap against the viscous flow of the gas: 
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(Lin & Papaloizou 1986). In the MMSN with a = 3 x 10~ 5 , M gap ~ 15M e at 5 AU. Numerical 
simulations show that gas flows onto a core even when a gap exists (D'Angelo et al. 2002), 
allowing a gas-giant planet to form. When M > M gap , the maximum envelope growth rate 
is 
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where Mj is the mass of Jupiter (Alibert et al. 2005; Veras & Armitage 2004). In the MMSN 
with a = 3 x 10~ 5 , a 10M e core at 5 AU will grow to Mj in ~ 4 Myr according to equation 7. 



3. Simulations of Growth with Migration 

We now examine the growth and migration of giant planets using a more detailed 
model. We consider a disk with 2.5 < a < 50 AU, containing 640 cores, each of 10~ 4 M e , 
and a population of planetesimals. Cores sweep up planetesimals following Eq. 1. At each 
radial location, e and i vary due to gas drag and perturbations from the cores, calculated 
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independently following Inaba et al. (2001) and Ohtsuki et al. (2002). Cores are assumed 
to have circular, coplanar orbits due to dynamical friction with planetesimals and tidal 
interaction with the gas. Neighbouring cores merge when their orbital separation < 8r^, 
maintaining the typical spacing of 10-12 seen in N-body simulations of oligarchic growth 
(Kokubo & Ida 1998). The capture cross section of a core's atmosphere is calculated following 
Inaba & Ikoma (2003). Cores undergo migration according to Eqs. 3 and 4. The maximum 
type-II migration rate is determined by the rate of angular momentum transport through 
the disk (Ida & Lin 2004). The core gap-opening mass is determined by Eqs. 5 and 6. 

The critical mass for a core to begin accreting a gaseous envelope is 



up planetesimals (Inaba et al. 2003; Ikoma et al. 2000; Ida & Lin 2004). Here, we assume 
k = 0.05 cm 2 g" 1 and M c = 20M e . When M > M crit , the envelope growth rate is 



(Ikoma et al. 2000; Ida & Lin 2004) where y = —1/4 is an empirical fit based on simulations 
by Hubickyj et al. (2005). For cores that have opened a gap, the maximum gas accretion rate 
is given by Eq. 7. Initially, £ so iid and S gas are given by Eqs. 2 with x = —1/2, modified by a 
factor of exp (—a/25 AU). This profile is shallower than the MMSN to allow for migration 
and planetesimal drift due to gas drag. The disk is gravitationally stable everywhere such 
that Q > 2. Initially S gas /S solid = 90, and gas disperses exponentially over time with a time 
constant of 1 Myr. Solid material has a density of 1.5 gem -3 , the stellar mass is 1M , and 
r and a are assumed to be independent of a and time. Simulations last for 10 million years. 

Figure 2 shows the evolution of three cores in a simulation with r — 100 m, a = 3x 10~ 5 , 
and S solid = 6 gem -2 at 5 AU. The upper panel shows the cores' orbital evolution, with arrows 
showing when each core opens a gap. The core at 5 AU grows rapidly and begins noticeable 
type-I migration after ~ 30, 000 years. At 50,000 years, the core opens a gap and type-II 
migration begins. However, inward migration is offset by the angular momentum gained by 
absorbing smaller bodies that migrate into the core's vicinity. These collisions can be seen 
as jumps in the core's mass, shown in the lower panel of Figure 2. After 40,000 years, the 
core starts to acquire an envelope. Envelope growth slows over time as E gas decreases and M 
increases, reducing accretion across the gap. A second core at 20 AU undergoes rapid type-I 
migration, then opens a gap when a ~ 9 AU. This core acquires an envelope but its growth 





where k is the opacity of the core's atmopshere, and M is the rate at which the core sweeps 
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lags behind the body at 5 AU. The core moves inwards by ~ 1 AU due to type-II migration, 
but migration slows when E gas becomes small. A third core at 23 AU grows and undergoes 
type-I migration, but remains too small to acquire an envelope before the gas disperses. 
These three planets are close analogues of Jupiter, Saturn and Uranus respectively. 

Figure 3 shows the outcome of three simulations with different disk masses. Each row 
of symbols shows the surviving planets, with symbol radius oc M 1//3 . The black and grey 
symbol segments show the solid and gas mass fractions respectively. The numbers to the left 
of the symbols indicate S so i i( j at 5 AU. The last row of symbols shows the giant planets of 
the Solar System, where solid mass fractions include elements heavier than helium in each 
planet's envelope. More massive disks lead to larger final planets and allow giant planets to 
form further from the star. The model predicts that gas-giant planets will form and survive 
migration provided that a<3x 1CT 4 and r < 0.5 km. Most simulations generate radially 
ordered systems: 1 or 2 gas-giant planets with 3 < a < 20 AU, followed by 1 or 2 large cores 
containing little gas, and finally a disk of sub-Earth-mass objects, akin to the Kuiper belt. 

Gas-giant planets produced in the model have solid-to-gas ratios similar to Saturn but 
higher than Jupiter. High solid fractions arise when the innermost planet absorbs other 
cores as they migrate inwards. In the model, closely spaced cores always coalesce. Jupiter's 
low solid-to-gas ratio suggests it gravitationally scattered nearby cores rather than absorbing 
them. In most simulations, several cores cross the disk's inner edge due to type-I migration. 
Presumably such objects would be lost into the star. These bodies contain only a few 
Earth masses and migrate rapidly, so their dynamical effect on nascent terrestrial planets is 
probably not severe. No planets are lost via type-II migration unless a > 10 -4 . A gas giant 
migrating through the inner Solar System would remove most of the solid material present. 
Terrestrial planets that formed subsequently would contain mostly ice-rich material from the 
outer disk (Raymond et al. 2006). The rocky compositions of the inner planets in the Solar 
System suggest that no giant planets were lost this way and that a was small. 

If a ~ 10~ 4 at the disk midplane, type-II migration times will be comparable to disk 
lifetimes. In a flared disk, the migration rate is roughly independent of a (see Eq. 4). Hence, 
giant planets should be common everywhere between their formation location and the inner 
edge of the disk. The model predicts that this is the case. Figure 4 shows the planets 
produced in 48 simulations with 100 < r < 500 m and 3 x 10~ 5 < a < 3 x 10 -4 . Giant 
planets are abundant at all distances from the inner edge of the disk out to 20 AU. This is 
in accord with the observed distribution of extrasolar planetary orbits (Marcy et al. 2005). 

I would like to thank Alan Boss, Lindsey Chambers, Fred Ciesla, Roman Rafikov and an 
anonymous referee. This work was supported by NASA's TPF foundation science program. 
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Fig. 1. — Growth of a core at 5 AU from the Sun according to Eq. 1, with e and % determined 
by an equilibrium between perturbations from nearby cores and damping due to gas drag. 
Migration is neglected. Here, S solid =10 gem -2 and S g as/S S oiid = 90. Each curve shows 
growth for a different r. The dotted line shows the type-I migration timescale. 
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Fig. 2. — Evolution of three surviving cores in a protoplanetary disk with r = 100 m, 
a = 3x 10~ 5 , and S so i id = 6 gem" 2 at 5 AU. Growth is calculated using the model described 
in Section 3. Upper panel: distance from the star. Arrows show when each core opens a gap 
in the disk. Lower panel: core masses. 
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Fig. 3. — The outcome of 3 simulations with r = 100 m and a = 3x 10~ 5 . The value of £ so iid 
at 5 AU in each case is indicated. Surviving planets are represented by symbols with radius 
oc M 1 / 3 . The black and grey segments of each symbol give the solid and gas mass fractions 
respectively. The symbols marked 'Sol Sys' shows the giant planets of the Solar System. 
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Fig. 4. — Planets generated in 48 simulations with 100 < r < 500 m and 3 x 10~ 5 < a < 
3 x 10- 4 . 
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